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The rupture of double-stranded DNA under stress is a key process in biophysics and nanotech¬ 
nology. In this article we consider the shear-induced rupture of short DNA duplexes, a system 
that has been given new importance by recently designed force sensors and nanotechnological de¬ 
vices. We argue that rupture must be understood as an activated process, where the duplex state 
is metastable and the strands will separate in a finite time that depends on the duplex length and 
the force applied. Thus, the critical shearing force required to rupture a duplex within a given 
experiment depends strongly on the time scale of observation. We use simple models of DNA to 
demonstrate that this approach naturally captures the experimentally observed dependence of the 
critical force on duplex length for a given observation time. In particular, the critical force is zero 
for the shortest duplexes, before rising sharply and then plateauing in the long length limit. The 
prevailing approach, based on identifying when the presence of each additional base pair within the 
duplex is thermodynamically unfavorable rather than allowing for metastability, does not predict 
a time-scale-dependent critical force and does not naturally incorporate a critical force of zero for 
the shortest duplexes. Additionally, motivated by a recently proposed force sensor, we investigate 
application of stress to a duplex in a mixed mode that interpolates between shearing and unzipping. 

As with pure shearing, the critical force depends on the time scale of observation; at a fixed time 
scale and duplex length, the critical force exhibits a sigmoidal dependence on the fraction of the 
duplex that is subject to shearing. 

PACS numbers: 87.15.A-, 87.15.-v, 87.14.G-, 87.14.gk 


I. INTRODUCTION 

Understanding the properties of stress-induced un¬ 
winding of a DNA duplex is crucial in analyzing many 
biophysical systems. Within the cell, enzymes exert large 
forces on DNA; for example in replication or transcrip¬ 
tion [1]. In the growing field of DNA nanotechnology, du¬ 
plexes within self-assembled structures and devices can 
be deliberately subjected to substantial stresses [2-6] to 
enhance the design possibilities. Recently, force-induced 
rupture of DNA has even been pioneered as a mech¬ 
anism for measuring forces [7-10]; for example, Wang 
and Ha used duplex rupture to estimate the force that 
a mechano-sensitive receptor in a cell experiences from 
a surface [8]. The importance of the response of DNA 
to stress and recent progress in single-molecule force- 
spectroscopy techniques [11, 12] have led to many de¬ 
tailed investigations of stress-induced duplex disruption 
[13-25]. Simultaneously, theoretical models have been 
developed to predict and explain the experimental results 
[13, 14, 17, 22-24, 26-36]. 

In this article we study the physics underlying shearing 
of short DNA duplexes, illustrated in Fig. 1 (a). We argue 
that the currently prevailing theory [23, 27, 28, 35, 36] 
fundamentally misinterprets the behavior of the critical 
rupture force as a function of duplex length. We first 
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frame our discussion with a toy model to highlight the 
basic physics, and review previous approaches to under¬ 
standing the problem. We then explore the system in 
more depth with a detailed coarse-grained DNA model 
(oxDNA). Our interpretation emphasizes the finite time 
scale of experimental observation, and predicts a crit¬ 
ical force that depends on this time scale as well as 
duplex length. OxDNA also allows us to explore the 
mechanosensor of Wang and Ha [8] , which involves a gen¬ 
eralization of shearing. 

A. Background theory: Basic models of shearing 

Shearing of a duplex involves applying antiparallel 
forces to either the opposite 3'-3' or 5'-5' ends of the 
bound duplex, as shown in Fig. 1(a). Such a setup has 
been explored experimentally [8, 13-15, 20-23] and theo¬ 
retically [13, 14, 22, 23, 27, 28, 32, 35, 36]. In early exper¬ 
iments [13-15, 20-22], loads were increased dynamically; 
more recently, Hatch et al were able to use a constant 
force [23]. The critical force f c above which the duplex 
ruptures is a key observable; Hatch et al. measured the 
shearing force for duplexes from length 12 to 50 bp (base 
pairs), and found that it increased with duplex length at 
short lengths before plateauing in the long length limit 
[23] (Fig. 2 (a)). Clearly, for any non-zero force, the stable 
state involves widely separated strands. At finite temper¬ 
ature, therefore, any duplex will rupture given sufficient 
time. Thus, the critical force should be understood as 
the force required to rupture a duplex within a given ex- 
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FIG. 1. Applying stress to short DNA molecules. Depending 
on where the two anti-parallel forces are applied, DNA can be 
subject to stress in different modes: (a) shearing mode, (b) 
unzipping mode and (c) mixed mode. In the shearing mode, 
the maximum extension per base pair gained from disrupting 
the duplex prior to complete strand separation (~ 0.3 nm) is 
much less than in the unzipping mode (~ 1.3nm), meaning 
that duplexes tend to unzip at lower forces than they shear 
[17, 23]. 


perimental time scale, and it is fundamentally important 
to view rupture under shearing in the context of time 
scales of experimental observation. For systems in which 
load is dynamically applied, the time scale over which 
the force changes is an additional complicating factor. 
We will restrict our discussion to the constant force case. 

At room temperature, opening of a duplex is widely 
believed to be a thermally activated process [30, 37-40]. 
The system must pass through an unfavorable transition 
state in which the strands are attached by a small number 
of base pairs; one would expect the rate to be exponen¬ 
tially suppressed by the free energy cost AG* of reaching 
this state from the fully formed duplex. Assuming a typ¬ 
ical free-energy cost of base-pair disruption of AGb p [41], 
AG* = TVAGbp — AGo where A is the duplex length and 
AGo is an TV-independent constant that incorporates, for 
example, the typical number of base pairs in a transition 
state. To a first approximation, the effect of a constant 
shearing force / is to modulate AG* as the transition 
state and duplex state have different extensibilities. 

rf 

AG*(/) = AAG bp -AG 0 + / (x A (f')-x t (f'))df , (1) 

in which x& (/) is the average extension of the fully- 
formed duplex under force /, and x t (f) is the extension 
of the transition state. Evaluating the integral in Eq. 1 
necessitates a detailed model of DNA, but physical in¬ 
sight can be obtained through the crude approximation 
of Jo (• x d(f ) - = -(NS - So)f , in which <5 is 

a positive constant representing the degree per base pair 
to which the disrupted duplex will tend to have a larger 
extension due to its single-stranded state. 5 q> like AGo, 
accounts for the number of base pairs in the transition 
state. 


Given the above description, a simple model for the 
rate 1 /r(f) of duplex rupture as a function of duplex 
length and force would be 

-F- = ko exp(—(AG*(/))/fcT) ^ 

= fc 0 exp(-(AAG bp - AGo - (NS - S 0 )f)/kT ), 

in which ko is an unknown rate constant. Thus for rup¬ 
ture to occur on a time scale of experimental observa- 
tion T obs , we require AG*(/)) < AGf. obs = fcTln(T obs £:o). 
This expression defines a critical force below which du¬ 
plexes will not tend to rupture during the experiment, 
and above which rupture is typical: 

_ AAG bp _ AG* obs + AG 0 
Jcy ’ NS-5 0 NS-6 0 ' 

At large A, / C (A) ~ AG^p/S and is A- independent. As 
A drops, the importance of the second term on the right- 
hand side grows and f c ( A) therefore decreases, reaching 
zero when AAGb p — AGo = AG* obs • Thus our simple 
model predicts that duplexes below a certain length will 
dissociate spontaneously on an experimental time scale; 
above this length the critical force required for rupture 
will rise from 0 and eventually plateau in the long length 
limit. Although the toy model is crude, the underlying 
physics is reasonable: the shearing stress acts to reduce 
the stability per base pair of the duplex state, allowing 
duplexes to rupture within a smaller observation time 
r 0 bs than at zero force. Long duplexes require a higher 
force than short ones, since if more base pairs are present 
then each individual base pair must be weakened further 
to result in the same barrier height AG* obg . At high 
enough forces fS > AGb p , breaking a base pair is ther¬ 
modynamically favorable and so even long duplexes will 
rupture rapidly, explaining the plateau in f c at large A. 

The behavior of Eq. 3 is illustrated in Fig. 2(a), in 
which the toy model is fitted to the experimental data 
of Hatch et ai [23]. For simplicity, we assume that the 
transition state is a single base-pair state; we thus set 
So = S and AGo = AGb p and fit AG^p/S and AG* obg /£ 
to the data. The observed behavior is clearly consistent 
with experiment. Further, the fitting parameters are rea¬ 
sonable given the crudeness of the model. If S ~ 0.15 nm, 
AGb p ~ 1.5 kcalmol -1 , similar to values reported by 
Ref. [41], and the barrier height AG* obg /AGb p ~ 6—7, 
implying that duplexes of 7-8 bp can dissociate within 
the observation time (Is), which is not unreasonable 
[38, 39, 55]. Importantly, as shown in Fig. 2 (b), the shape 
of f c ( A) depends on AG* obg . From Eq. 3, f c depends lin¬ 
early on AG* obs and hence (through Eq. 2) it decreases 
logarithmically with increasing r Q bs- This dependence is 
physically reasonable; for example, the maximal length 
of duplex that will spontaneously dissociate (rupture at 
/ = 0) in a time scale of milliseconds is clearly much less 
than the maximal length that will spontaneously disso¬ 
ciate in a time scale of years. It is also consistent with 
the fact that when stress is increased over time, rather 
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FIG. 2. Behavior of the critical shearing force, (a) Critical force f c as a function of duplex length as observed in experiment 
(circles [23]) and predicted by basic models of shearing. The toy model is shown as a solid line; a two parameter fit to the 
experimental data was performed, yielding A G* ohs /5 — 74.3 kcal mol -1 nm -1 and AGb P /£ = 10.9 kcal mol -1 nm -1 . The results 
are compared to those of deGennes’ original form [27] (dotted line) and the modified deGennes’ form [23] (dashed line), in both 
cases using the parameters obtained by Hatch et al. for fitting the modified form to the experimental data [23] . (b) Illustration 
of the dependence of the critical rupture force f c on AG* obs and hence the observation time as predicted by the toy model. 
AGj-Qbs, and hence r Q bs, increases from top to bottom (AG> obs G {4.4,10.4,14.4}kcalmol -1 ), and AGb P = 1.52kcalmol -1 , 
(c) Dependence of the rupture force on base pair stability AGb P , as predicted by the toy model. AGb P decreases from top to 
bottom (AGb P G {1.52,1.22,1.02}kcalmol -1 ), and AG} obs = 10.4kcalmol -1 . In (b) and (c), 6 = 0.14nm. 


than held constant, the duplex rupture force increases 
with the rate at which the force is applied [13-15, 20- 
22]. Nevertheless, regardless of how long the observation 
time, for infinite length chains the critical force will still 
plateau at / c (oo) = AGb p /£; this limiting behavior is 
simply reached more slowly with N for larger r 0 b s - 

In Fig. 2 (c) we show the effect of varying AGb p on f c . 
Reducing this base-pair stability leads to a lower plateau 
height f c ( oo) and sets the lower cutoff length below which 
f c = 0 to larger N. 

Broadly similar reasoning to that underlying the toy 
model has been applied to the unzipping of DNA (in 
which two strands are pulled apart from the same end of 
the duplex - see Fig. 1 (b)) [30]. Experiments in which 
force is applied dynamically have also been analyzed in 
terms of models of activated processes [13, 14, 22]. De¬ 
spite this history, the physical principles highlighted by 
the toy model have not been widely applied to under¬ 
stand the dependence of the critical shearing force on 
duplex length [23, 27, 28, 35, 36]. The alternative rea¬ 
soning originated from deGennes, who modelled DNA as 
a ladder with springs connecting both neighbors within 
a strand and bases paired by interstrand hydrogen bond¬ 
ing [27]. He calculated the mechanical equilibrium of 
this system under applied shearing stress, and posited 
that an individual base-pair spring with extension above 
a certain critical value would rupture. DeGennes showed 
that, for longer strands, the shearing stress can be spread 
out over several base pairs in mechanical equilibrium, re¬ 
ducing the strain at the duplex ends and allowing the 
duplex to withstand higher stresses than a single isolated 
base pair. His result for the critical force is [27] 

fc = 2/iX -1 tanh (xiV/2), (4) 

in which /i is the force required to rupture a single base 
pair and x = yZR/Q is a function of the spring con¬ 


stant between neighbors in a strand, Q, and the spring 
constant between base pairs R. Similarly to Eq. 3, Eq. 4 
also has an f c that increases with N and plateaus in 
the long length limit. The underlying physics is, how¬ 
ever, entirely different. The deGennes model assumes 
that rupture can only occur when each base pair is indi¬ 
vidually unstable; it does not allow for metastability and 
a finite time scale for rupture during which the system 
can break many base pairs to climb over a free-energy 
barrier. The shape of the curve is governed by x, a prop¬ 
erty of the duplex state, rather than the relative prop¬ 
erties of single-stranded and duplex DNA. A number of 
groups have explored extensions and improvements to the 
deGennes’ model [23, 27, 28, 35], but their approaches 
remain fundamentally based on the idea that the depen¬ 
dence of f c on N can be explained by the degree to which 
longer duplexes can spread mechanical stress over multi¬ 
ple base pairs, rather than because longer duplexes need 
a larger force to reduce the free-energy barrier opposing 
dissociation to a low enough value. 

As well as being based on fundamentally different 
physics, the two approaches also make qualitatively dis¬ 
tinct predictions. Most obviously, Eq. 4 predicts a finite 
rupture force for TV = 1, whereas Eq. 3 predicts f c > 0 
only for N > 7V 0 = (AG* obs + AGo)/AGb p . Further, 
Eq. 3 predicts that f c (N) depends logarithmically on the 
observation time r Q bs through its linear dependence on 
AGj-obs? by contrast, the deGennes theory does not in¬ 
clude r 0 bs- 

The experimental data from [23] in Fig. 2 (a) are sug¬ 
gestive of f c > 0 only above a certain value of N = 
TV 0 ~ 5 — 8 bp, as in fact was recognized by the authors. 
Nonetheless, they attempted to fit the deGennes model 
to their data [23], modifying Eq. 4 to 

fc = 2/iX _1 tanh (\(N - N open )/2 ) + 2/i (5) 
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which includes a new parameter 7V open and an additional 
term 2/i. The first parameter was set to N open = 7, in 
order to generate a finite Nq. They argued for a reduced 
effective duplex length N — 7 bp on the grounds that 
the end base pairs are always broken at room tempera¬ 
ture. However, such an approach is inconsistent with the 
thermodynamics of DNA as currently understood [41]; 
for example, hairpins with stems of three or four base 
pairs can be stable relative to the unfolded state at room 
temperature [41, 42]. Additionally, this approach does 
not include a dependence on r Q b s as we argued above; 
the curve is identical for experiments that take millisec¬ 
onds and experiments that take years. Finally, the value 
of x obtained from their fit, 0.147, is unphysically low. 
X = 0.147 would imply Q/R = 92.5, i.e. that stretching 
hydrogen bonds between base pairs is almost two orders 
of magnitude easier than extending the distance between 
stacked bases in a duplex. Such an unreasonable value, 
necessary to force the curve to reach large N values before 
plateauing, provides further evidence that the deGennes 
theory does not explain the shape of / C (7V). We note in 
passing that the reason for the term 2/i = 7.8 pN term in 
Eq. 5 as compared to Eq. 4 is unclear; this term shifts the 
curve slightly, but does not strongly influence the above 
discussion. 

Although the toy model is useful in providing under¬ 
standing, it is extremely crude and neglects a number 
of potentially important effects. Most significantly, tak- 
ing Jo (xd(f') ~ = -(NS - 6 0 )f in Eq. 1 is a 

very strong approximation. In this work we therefore use 
oxDNA [43-45], a nucleotide level coarse-grained model 
of DNA that captures both the elastic behavior of dou¬ 
ble and single strands, as well as the basic physics of 
bonding. It should therefore be able to describe in detail 
the effects of the shearing forces on a duplex, and cap¬ 
ture metastable states, which we will argue are critical in 
order to obtain an / C (7V) that is consistent with experi¬ 
ment. We note that oxDNA has the advantages of rep¬ 
resenting DNA structure and mechanics more accurately 
than simpler statistical models such as those of Refs. [27] 
and [28]. Further, oxDNA is not derived with either the 
assumptions of the toy model or of the deGennes’ model 
and its derivatives in mind. Rather, we use it to explore 
whether the findings of the toy model are robust to a 
more accurate model of DNA thermodynamics and me¬ 
chanics, and whether effects analogous to those identified 
by deGennes are likely to play a significant role. 


II. METHODS 

In oxDNA, each nucleotide is treated as a rigid body 
[43-45]. OxDNA nucleotides interact through potentials 
designed to mimic hydrogen-bonding, stacking, chain 
connectivity and excluded volume interactions; these in¬ 
teractions combine to allow the formation of right-handed 
double helices between complementary strands at low 
temperatures. OxDNA incorporates physically reason¬ 


able representations of the thermodynamics, mechanics 
and structure of single-stranded and duplex DNA [43, 44] 
— the key ingredients whose interplay is central to this 
system. The model has been shown to reproduce impor¬ 
tant aspects of basic processes such as hybridization [40] , 
toehold mediated strand displacement [46] and hairpin 
formation [47]. It has also been successfully applied to 
explore stress-induced transitions [48-50]. In this paper 
we use a sequence-averaged parameterization of oxDNA 
[43, 44], which is ideal for identifying generic trends. 

Since direct measurements of the rupture kinetics are 
very demanding, we estimate rates indirectly from the 
free-energy profiles of stressed duplexes. We simulate 
the model using the virtual move Monte Carlo (VMMC) 
algorithm of Whitelam and Geissler (the variant in the 
appendix of Ref. [51]). Shear stress is generated by 
applying two anti-parallel forces of magnitude / to the 
center of mass of the nucleotides at the 3' end of each 
strand (which are situated at opposite ends of the du¬ 
plex). We measure the free-energy profile as a function 
of the number of base pairs between the two strands (a 
base pair being defined by a hydrogen-bonding energy of 
less than 0.596 kcal mol -1 . A cut-off which is around 15% 
of the typical hydrogen-bonding energy in the model), us¬ 
ing umbrella sampling [52] to facilitate the measurement 
of less favorable states. We also use umbrella sampling 
to prohibit the strands from separating (once separated, 
the strands would never rebind), and for simplicity we 
only include hydrogen-bonding interactions between na¬ 
tive base pairs (those present in the fully-formed duplex). 
It was shown [40] that in oxDNA non-native base pairs 
have only a minor effect for non-repetitive sequences on 
the hybridization transition in the absence of force; the 
consequences of non-native base pairs for shearing in¬ 
duced rupture will be studied in detail elsewhere [53] . 

We obtain the free-energy profiles at different force val¬ 
ues, separated by 2 to 5 pN and spanning our range of in¬ 
terest, from umbrella sampling simulations at To = 23°C. 
The umbrella weights were adjusted iteratively to have 
a uniform sampling as a function of the number base 
pairs in the duplex. The starting configuration at each 
stage was chosen to be the final configuration from the 
previous iteration. At each force value, we ran 3 to 5 
independent simulations, each with approximately 10 9 
VMMC steps per particle. The resulting free energies 
from independent simulations agreed, confirming that 
our free-energy profiles are indeed converged. The value 
of the rupture force f c for a given AG| obs is then ob¬ 
tained via fitting the rupture force versus barrier height 
data to an interpolating piecewise cubic Hermite polyno¬ 
mial V; f c = V(A G| oba ). We determine the rupture force 
at other temperatures (T % To), from extrapolated free- 
energy profiles obtained using single histogram reweight¬ 
ing, based on the method of Ferrenberg and Swensden 
[54]. Note that the accuracy of the extrapolated free 
energies decreases when T is far from To. Therefore, 
we restricted our considerations to temperatures where 
|T — Tq| < 12 °C. 
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FIG. 3. Shearing and unzipping as represented by oxDNA. (a) Free-energy profiles of the bound duplex of length 21 base 
pairs at T — 23°C in the shearing mode. The pulling force increases from top to bottom: / shear £ {0, 30, 50, 70, 90}pN. (b) 
Rotation of the helix to maximize extension without disrupting structure. A 6-base-pair helix has its maximum extension when 
its helical axis is aligned with the stress; a 3-base-pair helix can rotate to increase its extension along the force, (c) Free-energy 
profiles of the bound duplex of length 21 base pairs at T = 23°C in the unzipping mode. The pulling force increases from top to 
bottom: / unzip £ {0,10,15, 20, 25}pN. (d) Rupture force as a function of duplex length in shearing and unzipping modes. The 
critical free-energy barrier height increases from top to bottom (AG* obs £ {3, 5,8,11} kcalmol -1 ). Larger AG* obg corresponds 
to longer experimental measuring times. 


III. RESULTS AND DISCUSSION 
A. Pure shearing and unzipping 

We measure free-energy profiles of shearing for a range 
of forces and duplex lengths; typical profiles are shown 
in Fig. 3 (a). At low forces, formation of each additional 
base pair is thermodynamically favorable but as the shear 
force is increased the slope of the profile becomes shal¬ 
lower and the duplex becomes less stable. On top of these 
overall trends, we see features that reflect the helical ge¬ 
ometry of the stressed duplex. A duplex that has intact 
base pairs equal to a half turn (z.e. approximately 6 to 7 
base pairs) has its maximum extension when the duplex 
axis is aligned with the force. When such a duplex loses 
or gains base pairs it can rotate its helical axis away from 
the direction of the force, allowing a greater extension per 
base pair (see Fig. 3 (b)). The result is a non-monotonic 
dependence of the free energy on the number of duplex 
base pairs; the effect is strongest for short duplex sec¬ 
tions, and explains the peak in free energy that appears 
at approximately 6 to 7 base pairs as the force increases. 
Such geometrical effects are absent in deGennes’ ladder 
model [27] and in most subsequent developments of that 
approach, although a similar effect in the context of a 2D 
ladder model was considered by Chakrabarti and Nelson 
[28]. 

For comparison, we also measure the profile for unzip¬ 
ping (when the two strands are pulled apart from the 
same end of the duplex, see Fig. 1(b)). Typical behav¬ 
ior as a function of / is shown in Fig. 3(c). As with 
shearing, base pairs become less stable as the force is in¬ 
creased. Two differences with shearing are clear. Firstly, 
the forces required to reduce stability are much lower; 
this is because a single unzipped base pair increases the 
extension of the DNA along the applied force much more 


than a single sheared base pair (this is evident in Fig. 1). 
Secondly, the non-monotonicity observed for shearing is 
absent for this geometry; regardless of the number of base 
pairs present it is always favorable to keep the helix axis 
perpendicular to the unzipping force. 

OxDNA clearly shows a more complex dependence of 
free energies on / than incorporated into the simple 
model. Nonetheless, we can identify a barrier height 
AG^m in each profile as the free-energy difference be¬ 
tween the (local) maximum in the profile with the small¬ 
est number of base pairs, and the lowest (local) mini¬ 
mum in the profile that has more base pairs than the 
maximum. We take these configurations to be proxies 
for the transition state and the duplex state, respec¬ 
tively, and assume that the dissociation rate is given 
by 1/t(/) oc exp(—(AGg im )/fcT). Previous studies with 
oxDNA in the absence of stress have shown that attach¬ 
ment rates of duplexes are only weakly dependent on 
overall duplex stability [40, 46]; this is consistent with 
the assumption that the detachment rate is largely de¬ 
termined by exp(—(AGg im )/fcT), as this barrier is the 
dominant factor in duplex stability. 

As with the toy model, we can define a rupture force 
f c that is the force required to reduce AGj im to a cer¬ 
tain value AG| obs , set by the time scale of experimental 
observation. We plot critical rupture forces for unzip¬ 
ping and shearing as a function of duplex length N in 
Fig. 3 (d) for a range of AG| obg values. Despite the ad¬ 
ditional complexity of oxDNA, the results are broadly 
consistent with the toy model; in particular, the depen¬ 
dence of / C (7V) on the observation time is similar. In Eq. 3 
the toy model predicts a linear dependence of / C (7V) on 
AG^bs (and hence a logarithmic dependence on r Q bs) at 
fixed TV; Fig. 4 (a) shows that oxDNA is in broad agree¬ 
ment. The agreement is not perfect, however, and in 
particular Eq. 3 predicts that the magnitude of the gra- 
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FIG. 4. The dependence of f c on AG* obs at fixed A as esti¬ 
mated by oxDNA. (a) Shearing mode, (b) Unzipping mode. 
The apparent non-monotonicity in the fitted gradient that is 
present for shearing but absent for unzipping is due to the 
number of base pairs in the transition state changing with 
applied force. 


dient of f c with AG| obg should decrease monotonically as 
A increases; this is not apparent in our fits as the gra¬ 
dient for A = 10 is shallower than for A = 13. In fact, 
this is due to non-linearities in the dependence of f c on 
AG| obs (when viewed at constant / c , the N = 10 curve 
is actually steeper). This non-linearity is inconsistent 
with the toy model, which assumes that the properties of 
the transition state do not change with force. It is clear 
from Fig. 3(a), however, the transition state moves to 
higher numbers of base pairs as the force increases and 
the geometrical effects discussed previously overwhelm 
the favorable base-pairing free energy for the final few 
base pairs. This can be interpreted as an increase in the 
offset parameter with /, which helps to explain why f c 
is not perfectly linear in AGj. obg . Unzipping lacks these 
kinds of geometrical effects; an equivalent graph of / C (A) 
against AG| obg in Fig. 4 (b) does not show this apparent 
non-monotonic dependency of the fitted gradient with A. 

The deGennes approach is based upon estimating the 
force at which each base pair is inherently unstable, so 
that rupture of the duplex occurs directly with each step 
being favorable. This criterion corresponds to the system 
ceasing to be metastable with a well-defined barrier. The 
toy model predicts that each base pair is inherently un¬ 



duplex length [bp] 



FIG. 5. The consequences of temperature and base-pair sta¬ 
bility. (a) Variation of / C (A) with temperature as predicted 
by oxDNA, defining / C (A) as the force required to reduce 
the free energy barrier, A G\ ohs /kT to 17 (corresponding to 
10 kcalmol -1 at 23 °C). (b) Fit of oxDNA model (squares) 
to the experimental data of Hatch et al. [23]. The fitting 
procedure generates a value for the rupture barrier height 
of AG* obs = 8.8 kcalmol -1 and includes a small rescaling of 
the temperature, for reasons described in the text. 


stable when f c = AGb p /£; this force is A-independent. 
Thus any variation in / C (A) with A as AG|. obg -A 0 can¬ 
not be explained by the toy model: it must arise from 
considerations such as those of deGennes, or geometric 
factors such as those highlighted above that are absent 
in the toy model. We will now argue that, for oxDNA at 
least, these factors only influence / C (A) at very low A. 

In Fig. 3(d) we show how / C (A) varies with AG^ obg . 
The lowest value of AGf. obg (corresponding to the short¬ 
est time scale) that we consider is 3kcalmol -1 ; smaller 
values make the unambiguous identification of maxima 
and minima in the free energy profile difficult, and in¬ 
deed the assumptions inherent in our kinetic model are 
less robust. Even at this finite value of AG* , how- 
ever, f c depends strongly on A only for A < 10; in the 
limit of AG| obg -A 0, this dependence will be even fur¬ 
ther truncated. Indeed, the shape of / C (A) at AG* obg = 
3 kcalmol -1 still reflects the predictions of the toy model 
and its metastability-based arguments (cutting through 
f c = 0 at finite A, for example). The relatively short 
range of A prior to the plateau as AG* obg -A 0 observed 
for oxDNA is consistent with the fact that geometric ef- 
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fects are only large for helices below a single pitch length, 
and that there is no physical reason why the x parameter 
in Eq. 4 should be much smaller than unity. 

OxDNA allows exploration of the effect of temperature 
changes. In Fig. 5 (a), we report / C (7V) at a range of tem¬ 
peratures. Both AG* and kT in the exponent of Eq. 1 
depend on T. Thus to compare critical forces at differ¬ 
ent temperatures for the same observation time, we find 
the force f c which gives AG\ irn /kT = 17 (lOkcalmol -1 
at 23°C, the standard temperature used in this study). 
We thereby include the full temperature dependence of 
the exponent in Eq. 2; any temperature dependence in 
the prefactor fco is neglected. Unsurprisingly, / C (7V) is 
reduced by temperature, as base pairs become less stable 
and free energy barriers of a certain height slightly easier 
to climb. Consequently the plateau force / c ( oo) drops 
and the length of duplex below which f c = 0, TVq, rises. 
These changes with temperature are similar to those seen 
for the toy model when base pair stability is varied, as 
shown in Fig. 2 (c) (note that in this case, the variation 
is performed at fixed AG* = AG* ob J. 

The data of Hatch et al were obtained at a lower salt 
concentration than those used to parameterize oxDNA 
to duplex formation in the absence of force 0.15 M 
as opposed to 0.5 M [Na + ]), implying weaker duplexes in 
experiment. Additionally, oxDNA is known to slightly 
overestimate (on the order of 10%) the critical force for 
the overstretching transition to ssDNA at 0.5 M [Na + ]; 
this is essentially the critical shearing force for infinitely 
long strands. As a result, the direct fit of oxDNA to 
the experimental data (using only AG* obg as a fitting 
variable) is only qualitative. To take these effects into 
account we fit to the data by shifting the temperature to 
T = 35°C, a slight (4% absolute change) increase when 
compared to the T = 23°C of the experiments. In ad¬ 
dition to the temperature, we also need to fit AG* obs 
which measures the effective barrier for rupture. Using 
these two fitting parameters, we find, as can be seen in 
Fig. 5 (b), very close agreement with the experiments of 
Hatch et al. The value of AG* = 8.8 kcalmol -1 is rea- 

1 obs 

sonable as it corresponds to approximately 6 or 7 base 
pairs [41], a sensible number for spontaneous rupture on 
the time scale of seconds [38, 39, 55]. Whilst the tempera¬ 
ture rescaling is crude, it is equivalent to a slight rescaling 
of the interaction strengths within oxDNA, producing a 
very similar model with slightly less stable base pairs. 
This oxDNA-like model, which still possesses mechanical 
and structural properties close to those of physical DNA, 
can quantitatively reproduce the experimental data. We 
therefore argue that the basic mechanism we propose, 
an activation based rupture of the duplex strands, can 
explain the experimental results. 


B. Mixed shearing and unzipping 

Given the success of oxDNA in describing the pure 
shearing mode shown in Fig. 1 (a) we next turn to the 


mixed shearing and unzipping mode shown in Fig. 1 (c). 
Wang and Ha have designed a tension gauge tether 
(TGT) to measure forces that involves applying stress 
in this “mixed mode”. In this setup, force is applied at 
the end of one strand and to any one base of the other 
strand, so that s base pairs are sheared and N — s + 1 
base pairs are unzipped. Note that the total appears to 
be N + 1: this is because there is no difference between 
unzipping and shearing of a single base pair; for the anal¬ 
ysis below it is helpful to define both quantities to include 
this ambiguous base pair. OxDNA allows us to explore 
the physics of this more complex system, using the same 
assumption that f c is the force required to reduce AG* obs 
to a specific value. Fig. 6 (a) shows the rupture force for 
a duplex of length N = 21 bp that is pulled in a mixed 
mode as a function of the number of sheared base pairs s. 
We see a characteristic sigmoidal curve for / < f nixed (A, s) 
at fixed N. f™ lxed (N, N) is identical to our pure shear¬ 
ing data in Fig. 3 (d) and depends on AG* c as expected, 
ymixed (TV, 1) is identical to the pure unzipping data. 

The sigmoidal curves therefore interpolates between 
the unzipping and shearing limits. A lower bound on 
the critical force is given by 

/ c mixed (V a) > max (/ c shear (s), / c unrip (i\0) , (6) 

as is clear from Fig. 6(b). The bound follows from 
two arguments. Firstly, the system can be viewed as 
a sheared duplex of s base pairs, with additional base 
pairs that must also be unzipped. Although these ad¬ 
ditional base pairs might not substantially stabilize the 
TGT at high force, they cannot act to destabilize it ei¬ 
ther; thus f™* ed (N, s) > / c shear (s). Secondly, sheared 
base pairs are always more stable than base pairs sub¬ 
ject to an unzipping force (Fig. 3(c)), due to the far 
greater increase in extension upon rupture in the sec¬ 
ond case. Therefore the system must be at least as sta¬ 
ble as an TV-base-pair duplex subject to pure unzipping; 

y mixed (^V, S ) > 

At large forces, the unzipped base pairs contribute al¬ 
most nothing to the barrier to dissociation and hence do 
not strongly affect the stability of the TGT, because at 
forces / > / c unzip (oo) each base pair subject to an un¬ 
zipping force is inherently unstable. Thus at large 8, 
ymixed (TV, s) follows J® hear (s) very closely as is evident 
in Fig. 6(b). On closer inspection, we note that in fact 
/™ixed(7V, 8) is slightly greater than /® hear (<s) for s < N. 
This discrepancy arises because in a purely sheared sys¬ 
tem, the base pairs present in the transition state can be 
anywhere within the duplex. In the TGT however, the 
unzipped base pairs at the end of the duplex anchor the 
final sheared base pairs (as illustrated schematically in 
Fig. 6(c)). This anchoring decreases the entropy of the 
transition state, and thus increases the barrier to TGT 
rupture slightly. 

Together, the lower bound and the insight that 
f™ ixed (N,s ) w / c shear (s) for / c mixed (iV,s) > / c unzip (oo) 
imply the sigmoidal curve for f™ lxed (N, s) at fixed N. 
At s = N, / c mixed (AT,s) = / c shear (Af). As s decreases, 
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FIG. 6. (a) Rupture force as a function of the number of sheared base pairs s , when a 21-base-pair duplex is pulled in the mixed 
mode. AG} obs G {5,8,11} kcalmol -1 , increasing from top to bottom, (b) Rupture force in the mixed mode as a function of 
sheared bases s for a duplex of length N = 21 (filled circles) (AG> obs = 10.4 kcalmol -1 ). (+) symbols show the corresponding 
rupture force for the pure shearing of duplex of length s and (x) symbols represent pure unzipping rupture force of a duplex 
of length N — s + 1. Open squares are Wang and Ha’s estimate of the rupture force [8]; their model is described in the text, 
(c) Schematic representation of the anchoring of the final sheared base pairs by the in the mixed mode by the base pairs that 
will eventually be unzipped; no such anchoring occurs for a pure shearing system. 


the curve f™* ed (N,s) approximately follows /c hear (s), 
which has the characteristic shape discussed previously 
and ym ixed (7V, s) drops towards 0 increasingly rapidly. 
Eventually, however, the lower bound f™* ed (N, s) > 
fu nzip(^y) comes i n t 0 play, forcing the curve to plateau 
at low s. 


Wang and Ha calibrated their TGT by assuming that 
/ c mixed (7V, s ) = / c shear (s) [8], with / c shear (s) given by Eq. 5 
(the formula of Hatch et al. [23]). They used the pa¬ 
rameterization of Hatch et al ., but set 7V 0 p en = 0 rather 
than TVopen = 7; due to this choice, /® hear (s) ^ 10 pN at 
5 = 1, rather than 0 (or 3.9 pN, as would be predicted 
by the original deGennes formula [27]). Their resultant 
calibration curve is shown in Fig. 6(b). Note that the 
shape is quite different from that predicted by oxDNA; 
in particular it has no inflection. The Wang and Ha curve 
predicts a physically reasonable f c for s ~ N by design, 
and, due to the choice of using the Hatch parameteri¬ 
zation of Eq. 5 but neglecting N open , predicts a roughly 
reasonable value for s = 1. In this picture, however, the 
finite f c at s = 1 is due to the robustness of a single 
sheared base pair rather than a duplex of length N that 
is subject to unzipping. Our results suggest that given 
agreement at s = 1 and s = N, the critical forces pre¬ 
dicted by Wang and Ha’s calibration approach would be 
far too high at intermediate values of s. Furthermore, 
we argue that f™* ed (N,s) depends on the time scale of 
observation; a consideration that does not appear in the 
Wang and Ha calibration curve. 


IV. CONCLUSION 

In this article we argue that shearing of short DNA du¬ 
plexes should be understood as an activated process, in 
which a metastable duplex dissociates into single strands. 
Whether or not a certain force is large enough to cause 
dissociation within an experiment therefore depends fun¬ 
damentally on the observation time r Q b s ; in principle, all 
duplexes subject to a constant force will dissociate per¬ 
manently given enough time. 

We used an extremely simple toy model to highlight 
the basic predictions of such an understanding. For a 
given r 0 bs 5 we must apply a large enough force to re¬ 
duce the free-energy barrier opposing dissociation, AG*, 
to AG* obg . Since more base pairs must be disrupted in 
longer duplexes to reach the transition state, the force- 
influenced stability per base pair must be less to achieve 
AG* = AG* obs ; hence the critical force f c increases with 
duplex length N. Once the force is large enough to re¬ 
duce the free-energy gain per base pair to zero, however, 
even the longest duplexes will have AG* < AG* obs . This 
introduces a maximum / c , leading to a plateau in / C (7V) 
for long duplexes. Importantly, this simple picture natu¬ 
rally predicts that / C (7V) tends to zero at finite N. This 
overall behavior is consistent with experimental observa¬ 
tions. 

To go beyond the simple model, we have used 
oxDNA, a coarse-grained model of DNA, to obtain more 
physically realistic force-dependent free-energy barriers. 
OxDNA reproduces the basic behavior predicted by the 
toy model, despite the more complicated free-energy 
profiles and force-extension properties. In particular, 
oxDNA also predicts that / C (7V) tends to zero at finite 
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TV, and that / C (7V) includes an approximately logarithmic 
dependence on the observation time. OxDNA is itself a 
simplified model, but the key aspects of DNA relevant to 
this study (force-extension properties, thermal stability, 
geometry) are known to be physically reasonable [43, 44]. 
The generic results are therefore likely to be robust; a 
detailed investigation into the relationship between rates 
and free-energy barriers for sheared oxDNA molecules, 
incorporating possible mis-aligned base pairs, is ongoing 
[53]. We have also demonstrated that, by slightly ad¬ 
justing the interaction strengths in oxDNA, in part to 
take into account the fact that it was parameterized at a 
different salt concentration, we can fit the available ex¬ 
perimental data quantitatively, supporting the essence of 
our claims. 

Previous work has contended that the shape of / c ( N) 
arises because a shorter duplex is immediately unsta¬ 
ble to rupture (rather than metastable with a short 
enough life time) at a lower force than longer duplexes 
[23, 27, 28, 35, 36]. We see evidence that shorter du¬ 
plexes are indeed unstable at lower forces in our data for 
oxDNA, but we argue that this is not the primary factor 
in determining the shape of / C (7V). Conceptually, it is 
clear that any theory that describes whether a molecular 
system undergoes change on a time scale of seconds or 
more should allow for metastability. In terms of concrete 
experimental evidence, it is known that duplexes of a few 
base pairs will spontaneously dissociate in the absence of 
force on a short time scale [38, 39, 55]. This is consistent 
with the data of Hatch et al. that show / C (7V) tending 
to zero at finite N [23]. This behavior does not natu¬ 
rally arise in models that predict / C (7V) based on absolute 
instability rather than metastability. Further evidence 
against the deGennes’ model arises from the observation 
that Hatch et al. needed to use unphysical values for the 
ratio of the spring constants between stacked neighbors 
in the same strand and hydrogen-bonded base pairs on 
opposite strands in order for the theory to fit their data. 

Further tests of the mechanisms that determine the 
critical force could be explored by experiments in which 


the observational time scale is varied. Hatch et al. did 
consider time scales between 1 s and 3 s, but given the 
approximately logarithmic variation expected in f c with 
t Qbs 5 it is unsurprising that systematic effects were not 
visible above experimental noise. The approach pre¬ 
sented in this work predicts a change in / C (7V) with r Q b s ; 
previous arguments based on absolute instability do not. 
Data from a wider range of duplex lengths N would help 
to differentiate between the typical shapes of / C (7V) pre¬ 
dicted by the two curves. We also make predictions for 
the temperature dependence of / c (TV), which could be 
tested in experiments. 

We have also used oxDNA to predict / ( f nixed (A/', s) for a 
mixed system in which s of the N base pairs are sheared, 
and the remainder are subject to an unzipping force. 
Such a system is the basis of the interesting tension gauge 
tether (TGT) proposed by Wang and Ha [8] for mea¬ 
suring biomolecular forces. We argue that the shape of 
ymixed(^y, g) j g ver y different from previously supposed, 
showing a characteristic sigmoidal shape as s is increased 
at fixed N. At small s, / < f nixed (A/', s) approaches the pure 
unzipping force /)) nzip (7V), but approximates the critical 
force for shearing an 5-base duplex, /c hear (s), at larger 
s. Not only do we predict that careful calibration of the 
TGT will reveal this complex / ( f nixed (A/', 5 ), we also claim 
that the time scale of experimental observation will in¬ 
fluence /™ lxed (iV, 5 ), in analogy with the pure shearing 
system. Thus quantitative use of the TGT will require 
extensive calibration. 
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